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The so-called $ parameter, which measures the transverse momentum or energy correlations 
(fluctuations) in high-energy collisions independently of the particle multiplicity, is computed for 
the equilibrium ideal gas. As expected, $ vanishes for the particles obeying Boltzmann statistics 
but is finite for the quantum particles, positive for bosons and negative for fermions. $ PJ _, which is 
Q\ ' found for the pions gas, significantly exceeds the value of <& p± measured by the NA49 experiment. 

t-H , The discrepancy is discussed. 
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There is a variety of correlations observed in proton-proton or proton-antiproton interactions at high energies. In 
particular, it has been found that the average particle transverse momentum depends on the particle multiplicity in 
a given collision . These correlations should be also present in nucleus-nucleus (A- A) collisions if such a collision 
is a superposition of nucleon-nucleon (N-N) interactions. However, there is no straightforward method to observe 
them since the final state particles in A-A collisions originate from the various N-N interactions while the correlated 
■ particles come only from the same N-N interaction. 

In our earlier paper || we have introduced a rather tricky quantity (later called <!>) which appears to be sensitive to 
the correlations independently of the particle multiplicity. If the A-A collision is a superposition of N-N interactions 
with no secondary collisions, the value of $ is exactly the same for the N-N and A-A case. If the secondary interactions 
play an important role in nucleus nucleus collisions, the correlation of interest is reduced Q and in the limiting case, 
when the final sate particles are independent from each other, <!> equals zero. The method developed in || has been 
recently applied to the experimental data and it has been found H| that the correlation, which is present in N N 
collisions, survives in proton-nucleus ones but is significantly reduced in the central Pb-Pb collisions. This result 
appears to be a very restrictive test of the models describing the N-N and A-A collisions. For example, the so-called 
random walk model Q is ruled out because it has been shown to produce, in contrast to the data, the stronger 
correlations in A-A than N-N case . 

Reduction of the correlations measured by $ in the central A-A collisions is naturally associated with the evolution 
of the system produced in these collisions towards the thermodynamical equilibrium. However, it has been correctly 
observed in Q that in the thermodynamical equilibrium the rudimentary correlation should be present. Our aim is 
to substantiate this observation. 

Let us first introduce the correlation (or fluctuation) measure <f> x , where x is a single particle characteristics such as 

the particle energy or transverse momentum. We define the variable z x = f x — x with the overline denoting averaging 
over a single particle inclusive distribution. As seen ~z x = 0. We now introduce the variable Z ', which is a multiparticle 

analog of z, defined as Z x = YliLi( x i — x )i where the summation runs over particles from a given event i.e. the 
particles which are produced in the collision. One observes that (Z x ) = 0, where (...) represents averaging over events 
(collisions). Finally, we define the measure $> x in the following way 
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Our purpose is to calculate $ x in the ideal quantum gas. At the beginning we identify x with the particle energy 
E and then we consider the particle transverse momentum p± . 
One immediately finds that 

-t_i r d*p.-. 2 _ i 



4 = - / ^(g-^ CTTi ■ ( 2 ) 



where the single particle average energy is 

E 

while p equals 



1 f d 3 p E 



p J (2tt) 3 \-ieP E ± 1 



9 J (2tt) 3 A-ie^ ± 1 ' 1 j 

/3 = T _1 is the inverse temperature; A = e^** denotes the fugacity and p, the chemical potential; E = m 2 + p 2 with 
m being the particle mass and p its momentum; the upper sign is for fermions while the lower one for bosons. It is 
worth noting that the result (Q) does not depend on the number of the particle internal degrees of freedom. 
Since Ze = U — NE, where U is the system energy, (Z E ) is computed as 
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S(V,T,A) , 



where H(V, T, A) is the grand canonical partition function || defined as 

S(y,T,A)=^^A Jv e-^ 



N a 



with V denoting the system volume and the index a numerating the system quantum states. As well known the 
grand canonical partition function of the quantum ideal gas equals 

d 3 p 



In E(y, T, A) = ±g V J j^L i n [ X ± X e~ 0E ] 



with <? being the number of the particle internal degrees of freedom. Consequently, 
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As previously the result is independent of g. One observes that 
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< z 2 and $ B < 
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for fermions, 



for bosons and 



(Z R > z | and $ s > 



(iV) 
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and $ B = 



in the classical limit i.e. when A -1 ^> 1. 

In the case of massless particles with vanishing chemical potential (which corresponds to A = 1), one finds $>e 
analytically. Namely, 
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and 
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where £(x) is the Riemann zeta function (C(3) = 1.202, £(5) = 1.037); the upper case is for fermions and the lower 
one for bosons. Further one computes 
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which finally give 



-0.07 
0.40 



T 



When the system is composed of the equal mass fermions and bosons with the numbers of the internal degrees of 
freedom gj and, respectively, g b , the analogs of the formulas (0) and (0) read 
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Using eqs. ^) and (Q, ||) one easily computes <J>_e for the baryonless quark-gluon plasma of two flavours where 
g f = 24 and g b = 16; <f> E = 0.17 T. 

Eqs. (|^, ^) can be used to get the measure $ Pi of the transverse momentum fluctuations. Since p± = psinO with 
p = |p| and 9 being the angle between p and the beam (z) axis one gets 
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where 
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Let us observe that $ p± is invariant under the Lorentz boosts along the beam axis. This is evident when eqs. (g, 
[Toj ) and (|[ |ll|) are written in the from which reveals the transformation properties. We consider as an example the 
average transverse momentum which can be expressed as 



P± = - 
P 
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with 



P = 



d 2 pj_dpn 1 

(2tt) 3 \-i e Pu" Pv ± i 



where u v is the four-velocity which determines the reference frame; u v = (1,0,0,0) corresponds to the thermostat 
rest frame. One sees that the two integrals, which determine p^, are both frame dependent due to the presence of 
dp\\. However, the dependence cancels out in the ratio of the integrals. Analogously one shows that $ Pi is invariant. 

We consider again the case of massless particles with vanishing chemical potential. As previously the calculations 
are performed in the thermostat rest frame. Then, 
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When the gas particles are massive and/or the chemical potential is finite, the correlation measure $ Pi can be 
numerically computed directly from eqs. (g) and ([l0|). In Figs. 1 and 2 we show $> p± as function of T and n for the 
pion gas. The pions are, of course, massive with m — 140 MeV. One sees that the presence of the finite mass reduces 
the correlation measure Q p± when compared to the massless case. 

As already mentioned Q p± has been experimentally measured in the central Pb-Pb collisions by the NA49 col- 
laboration. The result is <E> Pi = 0.7 ± 0.5 MeV [|[. If we identify the system freeze-out temperature with the slope 
parameter deduced from the pion transverse momentum distribution T = 180 MeV |hJ. Then, the value of & p±1 
which is read out from Fig. 1 for T — 180 MeV and \i = 0, equals 15 MeV and drastically exceeds the experimental 
value. The temperature is significantly reduced if the transverse hydrodynamic expansion is taken into account. The 
freeze-out temperature obtained by means of the simultaneous analysis of the single particle spectra and Bose-Einstein 
correlations is about 120 MeV (hJ. The value of $ Pi for T = 120 MeV and fi = equals 6.5 MeV and is still much 
larger than the experimental value. Let us discuss this puzzling discrepancy. 

& p± has been measured for pions which come form the limited phase-space region: 0.005 < pt < 1-5 GeV and 
4 < y < 5.5 §. However, it should not distort the value of $ Pi noticeably. First of all, one sees that the acceptance 
domain of pr covers the Pt— region which contributes to the integrals from eqs. (||) and (|Io|). One notes that the 
average pr approximately equals 2T. Secondly, we observe that the system longitudinal expansion influence the value 
of insignificantly as long as the transverse momentum distribution weakly depends on the particle rapidities. 
Finally, one notes that the size of the longitudinal momentum domain does not matter very much for the value of 
Q p± - Therefore, we conclude that the finite acceptance of the NA49 measurement cannot be responsible for the 
discussed discrepancy. 

We have considered three other ways to reconcile the experimental and theoretical values of . 

• The pions are out of chemical equilibrium and the chemical potential is negative. It appears however that /i 
must acquire an unrealistically large negative value (/j, = — 245 MeV) to get & p± =0.7 MeV. 

• A substantial fraction of the final state pions come from the hadron resonances. These pions do not 'feel' the 
Bose-Einstein statistics at freeze-out and should be treated as particles which obey the Boltzmann statistics. 
Then, they do not contribute to $ p± . Assuming that the fraction k of the final state pions come from the 
resonances, $ Pi is reduced approximately by a factor y/1 — k. The value of k must be again unrealistically large 
to reduce $ Pi sufficiently 
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• The Coulomb repulsion among the like-sign pions is known to significantly diminish the bosonic correlations, 
see e.g. However, taking into account the electromagnetic interaction should not change the value of $ p± 
noticeably. The point is that the Coulomb repulsion moderates the effect of boson statistics but the attraction 
among the unlike-sign pions generates the positive correlation of the similar size JTl| . 

It is possible that the combination of the three effects considered above sufficiently reduces the theoretical value of 
&p ± . However, it is not a simple problem to perform a numerically reliable analysis. Therefore, we leave it for the 
future studies. 

At the end let us comment on a somewhat paradoxical implication of our study. When the correlation measure <& 
was introduced ||, we expected that the value of <I> would be smaller in A-A than in p-p interactions. We shared a 
rather common opinion that the correlations observed in p-p case were of a dynamical origin and tended to be washed 
out in A-A by rescatterings. The recently obtained experimental data |^],[| comply with our expectation. Indeed, 
<5> Pi = 4.2 ± 0.5 MeV from p-p is about 6 times larger than that one from the central Pb-Pb collisions at the same 
energy 158 A • GeV However, we have now found that the equilibrium value of $> p± , which significantly exceeds 

$ p± from Pb-Pb, is close to that one from p-p. This implies that even in p-p collisions the origin of the correlations 
is poorly understood. Therefore, a complete solution of the puzzle raised by our equilibrium calculation requires a 
better understanding of the effects which control the value of <& at N-N level. 



I am indebted to Marek Gazdzicki for initiating this study and numerous fruitful discussions. The suggestion by 
Gunther Roland to consider the electromagnetic effects is also gratefully acknowleged. 
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Figure Captions 



Fig. 1. The correlation measure $ PJ _ as a function of temperature T for four values of the chemical potential fj,. The 
most upper line is for [i = 70 MeV, the second one for /i = etc. 

Fig. 2. The correlation measure $ Pi as a function of chemical potential /i for four values of the temperature T. The 
most upper line is for T — 200 MeV, the second one for T = 160 MeV etc. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9806089v3 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9806089v3 



